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We present methods of finding complex eigenvalues of the Liouvillian of an open quantum system. 
The goal is to find eigenvalues that cannot be predicted from the eigenvalues of the corresponding 
Hamiltonian. Our model is a T-type quantum dot with an infinitely long lead. We suggest the 
existence of the non-trivial eigenvalues of the Liouvillian in two ways: one way is to show that the 
original problem reduces to the problem of a two-particle Hamiltonian with a two-body interaction 
and the other way is to show that diagram expansion of the Green's function has correlation between 
the bra state and the ket state. We also introduce the integral equations equivalent to the original 
eigenvalue problem. 

PACS numbers: 



I. INTRODUCTION 

If a system has a compact potential in an infinite vol- 
ume, we can obtain complex eigenvalues of the Hamil- 
tonian under the boundary conditions of outgoing waves 
only or incoming waves only [H, , or with the Feshbach 
formalism 0, HJ- The eigenstates corresponding to the 
complex eigenvalues are not ghost states but have phys- 
ical meaning; they describe resonance states, which play 
crucial roles in high-energy physics, and the imaginary 
part of the complex eigenvalues is the inverse lifetime 
of the corresponding resonance state @, IH-12^. though 
they do not have probabilistic interpretation because of 
divergence of their norms (See, however, Ref. [1| for an 
extended probabilistic interpretation.). 

A discussion of the same sort can be applied to the 
Liouvillian, which is the generator of the time evolution 
of density matrices in the von Neumann equation: 



d 

i—p = Lp= [H,p] = zp, 



(1) 



where we use a unit with h = 1. In the infinite sys- 
tem with a compact potential, we expect that the Li- 
ouvillian can have complex eigenvalues z under certain 
boundary conditions, which have not been yet under- 
stood well, or with the Feshbach formalism for the Li- 
ouville operator. For finite systems, all the eigenvectors 
and the eigenvalues of the Liouvillian are written in the 



forms p n 



\n){m\ and z n 



E n - E m , respectively, 
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where \n) and E n are a discrete eigenstate and its eigen- 
value of the corresponding Hamiltonian, respectively For 
infinite systems, on the other hand, the discrete eigen- 
values of the Liouvillian may be of two types: one is 
z a /j = E a — En, where a and j3 are indices of discrete 
eigenvalues of the corresponding Hamiltonian including 
their complex eigenvalues, and the other is the type not 
written in the form E a — Ep, which we call nontrivial 
eigenvalues. An example of the nontrivial eigenvalue for 
the one-dimensional quantum gas interacting by the delta 
function potential has been discussed on the level of the 
weak coupling approximation by the one of the authors 
[23| . He showed that the nontrivial part of the eigenvalue 
gives a transport coefficient of the system. In contrast, 
we here present another example without any approxi- 
mation. 

In the present research, we present methods of finding 
nontrivial eigenvalues of the Liouvillian of an open quan- 
tum system. The reason why we are interested in such 
eigenvalues is that the eigenstates corresponding to them 
may include physically relevant non-equilibrium states, 
and the imaginary part of the eigenvalue can be the in- 
verse relaxation time. In contrast to finite systems, infi- 
nite systems have dissipation of particles into the infinity. 
Therefore, a state with relaxation can be an eigenstate of 
the Liouvillian under appropriate boundary conditions. 
Its eigenvalue may be complex and the corresponding 
eigenstate may be of the form ^ Q ^ c a a\a){P\ with some 
coefficients c a p. We can find such a state not by analyz- 
ing the Hamiltonian, only by analyzing the Liouvillian. 

The paper is organized as follows. In sec. [Hi we ex- 
plain a model that we treat in the present paper and 
our approach. In sees. IIIII and IIV1 we present methods of 
finding nontrivial eigenvalues of the Liouvillian and show 
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FIG. 1: The T-type quantum dot with an infinite lead. 

evidence for their existence. In sec. IIII1 we show it from 
the correspondence of the original Liouvillian and a two- 
particle Hamiltonian which has a two-body interaction. 
In sec. HVl on the other hand, we show it by indicating the 
operation of a Green's function on the bra state and the 
ket state simultaneously. In appendix \K\ we introduce 
some properties of the Liouville space and our notation. 
In appendix [B] we also show that the original problem 
reduces to two simultaneous integral equations. 

II. MODEL AND APPROACH 

Let us introduce the system that we consider. We an- 
alyze an open quantum system. The particle can then 
dissipate into the infinity. A T-type quantum dot with 
an infinitely long lead (Fig. [I} is the simplest model that 
has the above properties. The Hamiltonian is given by 
the following: 

oo 

H=-t (4+iCx + 4c+i) + Hi)(^+Cod)» ( 2 ) 

x — — oo 

where c\ and c x are the creation and annihilation oper- 
ators at the site x on the lead, while S and d are those 
at the dot site. We can rigorously obtain the spectrum 
of H . In particular, the Hamiltonian has complex eigen- 
values under the boundary conditions of out-going waves 
only or in-coming waves only. For this system, we define 
the projection operators P s and Q s as follows: 



P s = \d)(d\, 

OO 

Qs = ^2 \x)(x\, 



(3) 
(4) 



under appropriate boundary conditions for obtaining 
complex eigenvalues. We employ the Feshbach formalism 
0, H| , which was originally developed for the eigenvalue 
problem of Hamiltonians. Then, Eq. Q is written in the 
form 



L eS (Pp) = z(Pp), 



where 



L cS = PLP + PLQ- 



1 



-QLP 



z-QLQ 

with P and Q defined in Appendix |XJ We have 
PLP = (P s x P S ){H X 1 - 1 X H)(P S X P s ) 

= P S HP S X P s - P s X PsPPs = 



(10) 



(11) 



(12) 



because of Eq. (jSJ), where the operation 'x' is defined in 
Appendix [X] 

Since Eq. (TfOl) is given in the P space, we need the 
element ((d,d\L c g\d,d)): 



((d,d\L cS \d,d)} = tl((d,0\ 



1 



z-QLQ 



-tl((d,0\ 
-*?«0,d| 
+*?«0,d| 



1 



z-QLQ 



\d,0)) 
|0,d» 



1 



z-QLQ 



1 



z-QLQ 



\d,0)) 

|0,d», (13) 



where 



\d,d)) := \d)(d\ = tft|vac)<vac|(i, (14) 
0» := \d)(0\ = d f |vac)(vac|co, (15) 

and so on, with |vac) being the vacuum state. Thus, the 
problem of analyzing Eq. (jllj) is reduced to obtaining the 
Green's function of QLQ: 

QLQ = Q S HQ S x 1 - 1 x Q S HQ S 

+Q S HP S x Q s - Q s x P S HQ S 
+P s HQ s xQ s -Q s xQsHP s . (16) 



and thereby we have 

P S HP S = 0, 
P S HQ S = -hd^co, 
Q S HP S = -held, 



(5) 
(6) 
(7) 



Q s HQ s = -t (4+1^ + 4^+1)- (8) 

X — — oo 

Our aim is to solve the eigenvalue problem of the cor- 
responding Liouvillian: 



Lp = [H, p] = zp 



(9) 



III. CORRESPONDENCE OF THE 
ONE-PARTICLE LIOUVILLIAN AND THE 
TWO-PARTICLE HAMILTONIAN 

Let us introduce our main result of this paper in this 
section. We show that the one-particle Green's function 
of QLQ is equal to the two-particle Green's function of 
a new Hamiltonian %. In order to map the eigenvalue 
problem of the Liouvillian to the eigenvalue problem of 
the Hamiltonian, we use some powerful methods; i.e. the 
Dyson equation, closure, and so on. 

Though QLQ is a super-operator which acts on density 
operators, the density operators belong to a Hilbert space 



3 



with the inner product (IA4[) in Appendix [S] Therefore, 
we can treat QLQ as a usual operator which acts on the 
Hilbert space spanned by the states \i,j)). 

Let us here introduce a two-particle Hamiltonian: 

% = %lead + ^lead + %nt , (17) 

oo 

"ftfead = ~t (4 a *+l + a l + l a x), (18) 



The method of obtaining the one-particle Green's func- 



and (f\(Ei - Ul)- X \j) is 



tions (i'\(z — E± 
shown in Ref. [25(. 

We stress here that H has a two-body interaction Hi- 
This means that the evolution of the bra state and the ket 
state may correlate in the Liouville space, which suggests 
that some of the eigenvectors are not of the form |a)(/3| 
and their eigenvalues z are not of the type of E a — Ep . 



H b 



lead 



+t Oibw+i + bl +1 b w ), 



Hh 



-ti{a\a d + a d ao){l - b\b d ) 
-ti(blb d + b%)(l -a\a d ), 



(19) 



(20) 



where the a particle and the b particle are different and 
mutually commutative. Then the algebraic structure of 
H is exactly the same as that of QLQ in Eq. (JT5J): 



Q S HQ S x 1 
-1 x Q S HQ S 

Q S HP S x Q s 
-Q s x I'JIQ.. 

P S HQ S x Q s 
-Q s x Q S HP S 



u a 



lead 
6 

lead 

-tiala d (l - bifid) 
+tiblb d (l - a\a d ) 
-tia\a Q {l - bifid) 
+tib d b (l - a\a d ) 



In this way, the original problem of obtaining the Green's 
function of QLQ reduces to obtaining the Green's func- 
tion of the two-particle Hamiltonian H. 

Using the Dyson equation, we can rigorously show that 
the elements of the Green's function of H are represented 
by the Green's function of Ho, where we rearrange the 
Hamiltonian (|17[) as follows: 



H — Ho + Hi — H(j + Hq + Hi, 

oo 

HI = -t ^ (4 a x+i +al +1 a x ) 

X — — QO 

—ti(a]flQ + ala d ), 

oo 

ni = +t (bib x+ i+bi +1 b x ) 

X — — QO 

+ti(b% + blb d ), 
Hi = +ti(a' a d + a d ao)b d b d 
-ti(b f b d + b d b )a d a d , 



(21) 



(22) 



(23) 



(24) 



The two-particle Green's function of Ho is then given 
by the convolution of one-particle Green's functions as 
follows: 



(i'J\ 



1 



z - Ho 



\i,3) 



If 00 1 



- a M\ 



Ei - HI 



-|j)(25) 



IV. DIAGRAM EXPANSION OF THE GREEN'S 
FUNCTION 

In the present section, we show another sign of the 
existence of the nontrivial eigenvalues, which is not the 
type of E a — Ep , by describing the diagram expansion of 
the Green's function of QLQ. 

We now divide the Green's function of QLQ into the 
following two parts: 

QLQ=L + Li, (26) 

L =Q S HQ S x 1 - 1 x Q S HQ S , (27) 
Li =Q S HP S xQ s -Q s x P S HQ S 

+ P S HQ S xQ s -Q s x Q S HP S . (28) 

We then define G and Go as follows: 

1 



G 



Go 



z - QLQ' 
1 



z - Lo 



Then we obtain 



G — Gp y~^(XiGo) r 



(29) 
(30) 

(31) 



n=0 



using the resolvent expansion. 

Hereafter, we introduce our diagram exemplified in 

Fig.m 

• The upper line shows the time evolution of the bra 
state. 

• The lower line shows the time evolution of the ket 
state. 

• A thin vertical line indicates the action of Li, which 
move the particle from the lead to the dot or from 
the dot to the lead by one step. 

• A square indicates the action of Go- 

Considering that the necessary elements of G are those in 
Eq. (1131) , we can describe the term n — 2 of the expansion 
(f3"Tj) as Fig. (The diagrams expressed in the form of 
integral equations are given in Appendix [BT) 

The action of Go, which has information that the sys- 
tem has an infinite volume, acts on the bra state and 
the ket state simultaneously. This again suggests the 
existence of the nontrivial eigenvalues which cannot be 
described as E a — Ep . 
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FIG. 2: The diagrams for the term n = 2 of the expansion 

en. 



V. SUMMARY 

In the present paper, we considered the eigenvalue 
problem of the Liouvillian of a T-type quantum dot with 
an infinitely long lead and presented methods of finding 
nontrivial eigenvalues. 

In sees. IIIII and IIVI we showed evidence for the exis- 
tence of the non-trivial eigenvalues of the Liouvillian. In 
sec. IIIII we showed it from the correspondence of the orig- 
inal Liouvillian and the two-particle Hamiltonian which 
has a two-body interaction. In sec. HVI on the other hand, 
we showed it by indicating the simultaneous operation of 
Go on the bra state and ket state. Using the result in 
sec. IIIII we have obtained nontrivial eigenvalues approxi- 
mately by discretizing the integration in Eq. (|25[) . which 
may be reported elsewhere. 



Appendix A: The Liouville space 

In the present Appendix, we introduce some proper- 
ties and notations for simplicity of calculations. In the 
Liouville space, the operators are cumbersome to treat. 
The following notations are convenient for calculations. 

The Liouville space consists of density matrices p and is 
a Hilbert space with an inner product ({a' ft'\a/3)) defined 
as follows; 

\a t 0)) := \a)(0\, {\a)} : CONS, (Al) 
P = £3 C Q/3 1 £*,/?}), (A2) 

L\a,(3)) = [H,\a)(fi\] = (E a -E p )\a)(fi\ 

= (E a -E p )\a,p)), (A3) 

{{a',(3'\a,/3)) := tr[(|a')(/3'l) + (H</?|)] = <W <WA*) 

1 = ^\a,0))((a,P\. (A5) 

Throughout the paper, the operators that are denoted 
by capital letters without subscripts are in the Liouville 
space and operate on density matrices, while those with 



the subscript 's' operate on state vectors of the original 
Hilbert space. For simplification, we introduce the oper- 
ation ' x ' as follows: 

{A s x B s )p := A sP B s , (A6) 

A s x {bB s + cC s ) = bA s xB s + cA s x C s , (A7) 

(bB s + cC s )xA s = bB s xA s + cC s x A,, (A8) 

(A s xO)p = 0=(0xB s )p, (A9) 

{C s x D S )(A S x B s ) = C S A S x B S D S1 (A10) 

where p is a density matrix and b and c are c- numbers. 
Then the Liouvillian is given by 

L = H x 1 - 1 x H, (All) 

because then we have 

Lp = Hp- pH = [H,p\. (A12) 

We also define the projection operators 

P ■■= Ps x P s , (A13) 

Q := Q s xQ s + P s xQ s + Q s x P s , (A14) 

P + Q = Z, (A15) 

P s + Q S = I S , (A16) 

where operators denoted by I and 7 S are the identity op- 
erators in the respective spaces, P and P s are projections 
on the main system, and Q and Q s are projections on the 
environment. 

Appendix B: Integral equations 

In the present Appendix, we show that obtaining the 
elements of the Green's function of QLQ is equivalent to 
solving two simultaneous integral equations. 

Using the Dyson equation up to the second order, 



G — Go + GqLiGqLiG, 



(Bl) 



where the first-order term vanishes, we obtain the rela- 
tions: 

((x,d\G\0,d)) = ((x,d\G o \0,d)) 

-t\ J2((x,d\G \ Xl ,d)} 

Xl,X3 

(( Xl ,0\G \0,x 3 ))((d,x 3 \G\0,d)) 
+ tf J2(&d\Go\xi,d)) 

X\,X2 

({ Xl ,0\G \x 2 ,0)){{x 2 ,d\G\0,d)), (B2) 
«d >a ;|G|0,d»=t? ^2{{d,x\Go\d, Xl )) 

((0,x 1 \G Q \Q,x 3 ))((d,x 3 \G\0,d)} 
-t\ J2((d,x\G \d, Xl )} 

X\,X2 

{{Q, Xl \G \x2,Q)){{x2,d\G%d)), (B3) 
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({d,x\G\d,0}) = ((d,x\G o \d,0)) We define the Fourier elements of ((x, d\G\0, d)), 

4-/2 X^IM *\r Irf * \\ «d,a:|G|0,d», ((d, x\G\d,0)) , and ((a, d|G|d, 0)) as 

+ *i 2- {{*,x\<'o\d,x 1 )} G<M(k;d), G<°*>(d;k), G^ d0 \d;k) and G^(k; d) . Then 

3 we obtain two simultaneous integral equations: 
«0,a:i|Go|0,S3»«d,a:3|G|d,0» 

-t\ Y,{{d,x\G \d, Xl )) 

((0, nl Go 1 12, 0))((a; 2 , d\G\d, 0», (B4) 
<(M|G|d,0» = -tj ((x,d|G |xi,d)) 

«a 1) 0|Go|0,X3»«d,x 3 |G|d 1 0» 
+ t? ^((x,d|Go|.Ti,4 

Xl,X2 

«xi,0|Go|x2,0»((x 2 ,d|G|d,0». (B5) 



WOeQfi.. J\ _ 1 -t 2 r ^3, G {od) (d-k) , f 2 r (0d)(i„. j\ dk2 

u z +2tcosfe l lJ-i 2tt (z+2tcosfe)(z+2tcosfe-2tcosfc 2 ) 1 V ' ' J-tt 2tt (z+2t cos k)(z+2t cos fc-2t cos fc 2 ) 

G (d; k) — t-yG (d; fc) f_ w -^f ( z _2teosfc)(z+2t cosfci-2tcosfc) ~ *1 X-tt T^r" (z-2tcosfc)(z+2tcos fcl -2tcosfe) ' 



(B6) 



f((ffl)fj.i) _ 1 /2 dk 2 G (d0) (fc 2 ;rf) , ,2r<{dQi)(j. u\ f" dM 1 

^ V*:*} z-2tcosfc L 1J—!T 2tv (z-2t cos fc)(z-2t cos fc+2t cos fc 2 ) 1 V ' ^ > J —n 2ir (z — 2t cos fc)(z — 2i cos /c+2£ cos fc 2 ) ?T37\ 

WdOVi.. j\ _ 4.2 r (d0)( h . j\ f-w dfei 1 -i 2 r HS2... G (d0) (ri;fc 2 ) ^ ' 

i_z i-i "J" V^i^^J-tt 2tt (z+2tcosfc)(z-2tcosfci+2tcosfc) L lJ-n 2ir (z+2t cos k) (z~2t cos fc , +2t cos fc) - 



Then, the original problem reduces to solving the above solving the problem in the ways of sees. Mil and IIVI gives 
two simultaneous integral equations. In other words, a solution of the integral equations. 
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